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Abstract 

There is mounting experimental evidence that strong off-equihbrium phenom- 
ena, such as "memory" or "aging" effects, play a crucial role in the physics of 
vortices in type II superconductors. In the framework of a recently introduced 
schematic vortex model, we describe the out of equilibrium properties of vortex 
matter. We develop a unified description of "memory" phenomena in magnetic 
and transport properties, such as magnetisation loops and their "anomalous" 2nd 
peak, logarithmic creep, "anomalous" finite creep rate for T — > 0, "memory" and 
"irreversibility" of I-V characteristics, time dependent critical currents, "rejuve- 
nation" and "aging" of the system response. 
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1 Introduction 



The properties of vortex dynamics in type II superconductors crucially affect the 
overall system behaviour and have, thus, relevant effects in technological appli- 
cations Ijl], 0, H, H . In particular, in the last few years it has been discovered that 
vortex matter exhibits important, even dominant, history dependent phenomena 
in magnetic and transport properties, such as memory and hysteresis in mag- 
netisation curves along with irreversibility and aging in I-V characteristics (see 
[|, |, |3|, H, g5H3|, H, H, gl], 13, H, H, 0, |6| and references therein). 



These phenomena are markedly out of equilibrium effects and, here, we discuss 
their features as they emerge from the off equilibrium dynamics of vortex matter 
(see Ref.s in i |, 0, |, |). 

The above experimental findings have interesting analogies with "memory" 
and "aging" effects observed in other glass formers, such as polymers, supercooled 
liquids or random magnets P, H. Interestingly, "glassy" dynamics have important 
universal structural properties P, 0. Off equilibrium features arise when typical 
experimental probing times get much shorter compared to the system long (often 
inaccessibly long) intrinsic relaxation time scales. These can become huge at low 
temperatures or high densities, where a true equilibrium glass phase transition 
in some cases can be also found. Such glass transitions are called "ideal" 
because, as just stated, equilibrium might be hardly approached. In facts, the 
notion of "glassy phases" has been repeatedly used in relation to new equilibrium 
phases of vortex matter |10|, |ll], |13|- We are concerned here, however, 
with the general properties of off equilibrium dynamics of vortices, not with their 



equilibrium transitions [14 



We consider a schematic model [|T5] that contains the essential degrees of free- 
dom of a vortex system and is simple enough to allow a complete understanding of 
its off equilibrium dynamics in the same perspective successfully used for other 
glassy systems |^. The model (a coarse grained [p!^, p!7[] system of repulsive 
particles wandering in a pinning landscape in presence of a thermal bath and 
an external drive) describes several phenomena of vortex physics, ranging from a 
reentrant phase diagram in the {B, T) (field-temperature) plane, to the anoma- 
lous "second peak" in magnetisation loops (the "fishtail"), logarithmic creep and 
"aging" of magnetic relaxation, the finite creep rates for T — »• (without use of 
"quantum effects"), "memory" and history dependent behaviours in vortex flow 
and in I-V characteristics, and many others [T^ . 

We describe here the properties of such a model and depict a unified picture 
of creep and transport measurements. In particular, the system dynamics can 
be described by identifying its important time scales and their dependence on 
temperature, magnetic field and applied electrical current. We also suggest new 
experiments that will help to clarify the nature of glassy aspects in superconduc- 
tors. 

In the next section we introduce the model [|1^ and in the following we sys- 
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tematically compare its behaviours with experiments on magnetic and transport 
properties. Finally in the conclusions we give an overview of our scenario of off 
equilibrium phenomena in vortex matter. 



2 The R.O.M. Model 

Vortices in type II superconductors are described by the Ginzburg-Landau equa- 
tions. The typical high vortex densities and long interaction range imply that 
the vortex system is strongly interacting. In brief, this makes the theoretical 
description of its equilibrium and, even worse, dynamical properties highly non 
trivial 

An appealing and much used approximation for the microscopic vortex dy- 
namics is based on Molecular Dynamics (MD) simulations (see for instance Ref.s 
in [|l^, |19|, |2l|). However, even this simplified approach is hardly feasible to 
explore the physics of the long time and space scales, low temperatures and high 



densities region where glassy features substantially appear Alternatively 
it was proposed to use schematic discrete time and space models [^] to study 
vortex properties. 

More generally, to describe the relevant degrees of freedom of the vortex sys- 
tem one can introduce useful coarse graining methods, successfully applied to deal 
with many other multiscale problems (such as magnetism or crystals defects, see 
Ref.s in [^]). For clarity, let's consider the simple case of a system of straight 



parallel vortex lines, corresponding to a magnetic field B along the 2;-axis, where 
vortices interact through a two-body potential 0: 

A{r) = ^ [Ko(r/A') - K,{rlO] , (1) 

Kq being the MacDonald function, ^ and A the correlation and penetration 
lengths (^' = c^/a/2. A' = cA, c = (1 — B / Bc2)~^^'^)- A simple apphcation of 
the above methods in the present case, proposed in [jl^, |l5l, consists in coarse 
graining the vortex system in the xy-plane by introducing a square grid of lattice 
spacing, /q, of the order of the London length, A (see Fig.|I|). 

By this procedure, the original vortex system is mapped into a lattice model 
characterised by a classical field, n^, representing the number of vortices on the 
i-th coarse grained cell (see Fig.|I]). The presence in superconductors of an upper 
critical field, i?c2, implies that Ui must be an integer number smaller than Nc2 = 



Bc2lo/4>o nil (00 = hc/2e is the flux quantum). The Hamiltonian of the coarse 



grained model is [p!5|: 



2^ ' ' 2 



The first two terms of 7i describe the repulsion between the vortices and their 
self energy, and the last the interaction with a random pinning background. 
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Figure 1: A schematic plot of the procedure introduced to define the ROM lattice 
model. The original vortex system (left), coarse grained in "cells" of size Iq, is 
mapped into a lattice field model (right). 



For sake of simplicity, since Iq ~ A, we can consider the simplest version of 
Ti: we choose An = Aq = 1; Aij = Ai < Aq if i and j are nearest neighbours; 
Aij = otherwise; the random pinning is taken to be delta-distributed P^A^) = 
(1 — p)6{A'^) + pS{A''' — Aq) (see [^). We express all energy scales in units of 



Aq and, in particular, consider the important ratio k* = Ai/Aq. The existence 
of two possible orientations of the vortices can be taken into account by giving 
the particles, rij, a "charge" Sj = ±1 |T|, 0. Neighbouring particles with opposite 
"charge" annihilate. The external applied field controls the overall system "charge 
density" and thus a chemical potential term — /i X]j SiUi must be added to the 
Hamiltonian in Eq.(0) (where rij is replaced by s^nj). 

A standard mean field replica theory [1^ allows to evaluate the equilibrium 
phase diagram in the field-temperature plane of the above Hamiltonian, as shown 
in Fig.^ In absence of disorder it has, at low temperatures, a reentrant order- 
disorder transition in agreement with predictions and experiments on vortices 
in superconductors (see Ref.[jl], |], |] or, for instance, data on 2H-NbSe2 super- 
conductors from Ref. ||32||). For moderate values of the pinning energy {Aq < Ai), 
a second order transition still takes place, which at sufficiently strong pinning 
is expected to become a "glassy" transition, as is seen in Random Field Ising 
Models [0. The extension of the low T phase shrinks by increasing Aq (i.e., the 
highest critical temperature, T^, decreases) and the higher is k* the smaller the 
reentrant region (facts in agreement with experiments, see for instance Ref.s in 

^, 0). The above phase diagram can help to compare experimental and 



model temperature/field scales. 

We now go beyond mean field theory and discuss the dynamics of the model. 
First we consider the case where the external current is absent, i.e., there is 
no Lorentz drive on vortices. The simplest consistent approach to simulate the 
system relaxation at non-zero temperatures is a Monte Carlo Kawasaki dynamics 
on a square lattice of size L at a temperature T (see |^). This is a very 
standard approach in computer simulations of dynamical processes in complex 
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Figure 2: Main frame The mean field phase diagram of the ROM model in the 
plane {H*,T*) (T* = T/Ai and H* = fi/ksT are the dimensionless temperature 
and chemical potential of the applied field), for k.* = 10 and Aq = 0.0; 0.5; 0.75 
(res. full, dotted and dashed lines) and k* — 3.3 and = 0.0 (long dashed line). 
Inset The magnetisation Bean profile, M{x), as a function of the transversal 
spatial coordinate x/L (L is the system linear size), (for the shown values Next), 
in the 2D ROM model (T = 0.3, 7 = 1.1 lO^^). Notice the change in shapes 
when Next crosses Ngp — 13.5 (filled v.s. empty symbols). 
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fluids [^. In particular, we consider a system periodic in the y-direction. Its 
two edges parallel to the y-axis are in contact with a vortex reservoir, i.e., an 
external magnetic field, of density Next- Particles can enter and leave the system 
only through the reservoir. 

The above model, called ROM (Restricted Occupancy Model), is described 
in full details in |jl5[. It is extremely schematic, thus, also fully tractable, and, 
interestingly, it is able to describe many of the experimental observations on 
magnetic and transport properties of vortex physics. 



K*=0.26 7=2.2 10" 
O K*=0.26 7=1.1 10'^ 
^"=0.10 7=1.1 10"^ 




ext 



Figure 3: Main frame The magnetisation, M, as a function of the applied field 
density, N^xt, in the ROM model at T = 0.3 for the shown sweep rates 7 and 
K*. Notice the appearance of a "second magnetisation peak" when k* is large 
enough. Inset The equilibrium value of M (i.e., when the field ramp rate 7^0) 
at T = 0.3 {k* = 0.26). 



3 The Magnetisation 

The simplest quantity to characterise the vortices system is the magnetisation, 
which we now consider. In this section we will draw a picture of time dependent 
magnetic features such as magnetisation loops, their 2nd peak, "aging" creep, as 
well as phenomena like the finite creep rate, Sa > 0, found when T —>■ 0. 

The system is prepared by zero field cooling at a given T and then increasing 
the external field. Next, with a constant rate, 7. During the ramp of Next we 
record the magnetisation 

M{t) = Nin{t) - Nextit) ■ (3) 

Here A^j„ = and the Monte Carlo time, t, is measured in units of com- 

plete Monte Carlo lattice sweeps. 
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3.1 Magnetisation loops 

At low temperatures pronounced hysteretic magnetisation loops are seen when 
M is plotted as a function of N^xt (see Fig. Furthermore, when the parameter 
K* = Ai/Aq {k* can be directly related to the Ginzburg-Landau parameter k = 
X/^ mH) is high enough, a definite second peak ("fish-tail") appears in M. Very 
similar magnetisation data are observed in a number of different superconductors 
from intermediate to high k values (see, for instance, 0, ||, 0, ^ |3^, |35|, ^ 
, ^ and references therein). 
The actual shape of loops depends on the system parameters (and its size). 
In particular, the sweep rate of the external field, 7, is very important, as shown 
in Fig.|^. As soon as the inverse of the sweep rate is smaller than the system 
characteristic relaxation time (see below) strong hysteresis effects are present. 
Although the second peak does depend on dynamics through 7, in the ROM 
model it is related to a new phase transition: in the 7^0 limit, its location, 
Nsp, is associated with a sharp jump in Meg = lim^^o ^(7) (see inset of Fig.^). 
These findings are consistent with experiments (for instance, see Ref . |]32|, |33| , ^ 
P7| , P^ , |iD[1 ) and to some extents reconcile previously proposed 

opposite descriptions ("static" v.s. "dynamic"). 

3.2 History dependent relaxation: "aging" creep 

The presence, at low temperature, of sweep rate dependent hysteretic cycles, 
slowly relaxing magnetisation, and similar effects, indicate that our system, on 
the observed time scales, is not at equilibrium. We turn, therefore, to the theoret- 
ical description of the system dynamics by investigation of two times correlation 
functions. At the given working value of the applied field, we also record the mag- 
netic correlation function, C(t,t^) (with t > t^), which gives richer information 
than M{t) ||: 

C(t,t^) = ([M(t)-M(t^)]2) . (4) 

At not too low temperatures, for instance at T = 1.0 (a comparison of such 
a T value with experimental scales can be derived from Fig. |^), the system 
creep is characterised by finite relaxation times and no "aging" is seen: C(t,t^) 
is a function of t — t^. At long times, Cit.t^) is well fitted by the so called 
Kohlrausch-Wilhams- Watts (KWW) law p|: 

C(t,U^Coo{l-e-[(*-*-)/-]'} . (5) 

Eq.(|^) defines the characteristic time scale of magnetic relaxation, r. This is a 
crucial quantity to be considered when dealing with dynamical aspects of mag- 
netisation. The Kohlrausch-exponent, /?, and r strongly depend on T (a fact 
to be discussed below, see Fig.||) and on the applied field N^xt [|15|. The pre- 
asymptotic dynamics (i.e., t — << r) is also interesting and characterised by 
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Figure 4: Main frame The system equilibration time, r, from eq.(|^), enor- 
mously grows by decreasing the temperature T (here N^xt = 10). Below the 
crossover temperature Tg ~ 0.25, r is larger than the observation window. Inset 
Close to Tg, r plotted as a function of 1/T approximately shows a Vogel-Tamman- 
Fulcher behaviour, the continuous line (see eq.(^). 

various regimes. In particular, for not too short times, a power law is observed 
over several decades. 

The scenario described for T = 1.0 is found in a broad region at low temper- 
atures. However, around T = 0.5, a steep increase of r is found (see Fig.|). For 
instance at Next = 10, for temperatures below Tg ~ 0.25, the characteristic time 
gets larger than our recording window and the system definitely loses contact 
with equilibrium. The crossover temperature, Tg{Next) (which may be a func- 
tion of 7) has a physical meaning similar to the so called phenomenological glass 
transition point in supercooled liquids |^]. The presence of an underlying "ideal" 
glass transition point, Tc{Next), is often located by some fit of the high T data 
for r (see inset of Fig.^), such as a Vogel-Tamman-Fulcher (VTF) (or a power) 
law: 

- = -oexp(^). (6) 

Our data in 2D, are consistent with = 0. Interestingly, the VTF behaviour 
found here is in agreement with results from Molecular Dynamics simulations 
of more realistic London-Langevin models [|T^, [TB|, |2D|, In particular, 

the analogies with "window glasses" have been also outlined in the first of Refs 

. A VTF behaviour has been already experimentally observed in measures on 
samples resistivity (see [p^ ). 

Since below Tg relaxation times are huge, one might expect that the motion 
of the particles essentially stops. Instead, as shown below, the off equilibrium 
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Figure 5: Inset Logarithmic time relaxation of the two-times correlation function, 
C(t, tw), as a function of t — for the shown values of ty,. Data are recorded at 
T = 0.1 (i.e., below Tg) and Next = 16. The continuous lines are logarithmic fits 
from eq.(l^). Main Frame Off equilibrium dynamical scaling. The relaxation 
data of C{t, t^) from the inset and those recorded at N^xt = 4, 10, 16 for each of 
the shown are superimposed on the same master function. The asymptotic 
scaling is C{t,tw) ~ Sit/t^). 



dynamics has remarkably rich "aging" properties. In the inset of Fig.|^ we show 
that C{t,tuj), at T = 0.1, exhibits strong "aging": C depends on both times t 
and t^; in particular, its evolution (see Fig.|) is slower the older is the "age" 
("stiffening"). In the entire low T region (T < Tg), after a short initial power law 
behaviour, C{t, t^) can be well fitted by a generalisation of a known interpolation 
formula, often experimentally used [jl], |^, which now depends on the waiting time, 



C(t, ty 



1 + 



In 



tw + ^0 



(7) 



Eq. (^ (in agreement with the general scenario of Ref. |^5| ) implies the presence 
of scaling properties of purely dynamical origin (see Fig.^: for times large enough 
(but smaller than the equilibration time), C is a universal function of the ratio 

t/tyj-. C{t,tyj) S{t/tyj). 

In experiments about vortex creep in superconductors a crossover is usually 
found from a low T region with logarithmic creep to an high T region with 
typically power law or stretched exponential relaxations (see for instance [0 and 
Ref.s in [^). In particular, aging in magnetic creep has been recently observed in 
BSCCO samples |^6[. We also recall that the above phenomena are intriguingly 
common to many different systems ranging from polymers, to supercooled liquids 
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Figure 6: The vortex mean square displacement -R^(t) at Nf.xt = 10 for several 
temperatures. Below Tg ~ 0.25, B^{t) is strongly sub diffusive: B^{t) ~ f with 
z/ < 1. Straight hues are guides for the eye. 



0, spin glasses granular media p3| . 



3.3 Vortex mean square displacement 

The microscopic origin of the above features in the system dynamics can be 
understood by considering the vortex mean square displacement, R^{t) (plotted 
in Fig.|]for N(.xt = 10). At high enough T, B?{t) is linear in t (in agreement with 



experiments and MD simulations, see |^T], ^ and ref.s therein), but at lower 
temperatures it shows a pronounced bending. Finally, below Tg, the process 
becomes strongly sub diffusive: 

R^{t) ~ f (8) 

with V « 1. From this point of view, Tg is the location of a sort of structural 
arrest of the system, where particles displacement is dramatically suppressed. 
Each vortex is caged by other neighbouring vortices for long times. The system 
dynamics needs large scale "cooperative rearrangements" to relax [^. Interest- 
ingly, a very similar scenario has been recorded in real superconducting samples 



fsee for instance [4S 



3.4 The creep rate, Sa, for T ^ 



With the insight on the system dynamics obtained in the previous sections, we 



can now understand an other intriguing experimental observation ||53, |5^, p6 



about vortex matter: even at very low temperatures (where activated processes 
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Figure 7: Main frame The creep rate, Sa, in the ROM model for N^xt = 10 
as a function of the temperature, T, in units of Aq {k* = 0.28, 7 = 10~^). The 
superimposed hne is a hnear fit. Inset Creep rate, S, in a BSCCO single crystal 
at 880 Oe (from Aupke et al. [53]). 



should be absent) magnetic relaxation does take place. This surprising phe- 
nomenon, previously interpreted in terms of "quantum tunnelling" of vortices 
0, is also found in the present purely "classical" vortex model. More generally, 
we show here that a non-zero creep rate for T — is to be expected in systems 
"aging" in their off equilibrium dynamics. 

Experiments investigate the temperature dependence of the creep rate, tS^j, 
(see Fig. |^, where 

91n(M) 



Sa 



ain(t) 



(9) 



{Sa is, as usual, averaged in some given temporal window p, ^ |56|). 

When the temperature is extremely low the magnetisation still logarithmically 
relaxes (see inset Fig.^), and in both experiments and in our simulations, Sa 
approaches a finite plateau, SaiS^) > 0, for T ^ 0. In Fig.^ we plot the creep 
rate, 5*^, as a function of T. For comparison we present experimental data in 
BSCCO (from Ref.|5^) as inset (note that the values of Sa in our model and in 
real samples are very similar). In particular, we find that a linear fit of Sa{T) in 
the low T regime is very satisfactory (see Fig.|^): 



Sa{T) = 5° + aT 



(10) 



where both 5*^ and a are functions of the applied field Next- In the present model, 
as much as in experiments |^ |^, Sa(T) is non monotonous in T: at high T it 
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starts decreasing (this is due to the fact that, for a given observation window, at 
higher T the system gets closer to equihbrium, see Fig.^ and 
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Figure 8: Main frame The T ^ hmit of the creep rate, Sa, in the ROM 
model as a function of the applied field N^xt (for T = 10~^ and 7 = 10~^). 
The superimposed dashed curve is a power law to guide the eye. Inset: The 
relaxation of the magnetisation, M{t), in the model for Next = 10 and T = 0.25 
as a function of time (7 = 10^'^). The continuous line is the logarithmic fit of the 
text. 

By varying the applied field we find a range of values for 5"° very similar to 
experimental ones |5^, ^ |56l (see Fig.H). In particular, seems to decrease 
on average by increasing the field Next- The overall behaviour can be roughly 
interpolated with a power law: S'^{Next) — (Next/No)'^, where, for k* = 0.28, 
A^o — 0.01 and x ~ 0.6. As shown in Fig.|^, the presence of a small exponent, 
X, implies that sensible variations in Sq can be seen only by changing Next of 
orders of magnitude. Note that in Fig. || the dips in the Sa{0) versus Next data 
found at certain values of Next (namely around 3, 13, and 20) are statistically 
significant. They are located respectively close to the region of the low field 
order-disorder transition (see Fig.^, the 2nd peak transition and the reentrant 
high field order-disorder transition. 

In the slow off equilibrium relaxation at very low temperatures no activation 
over barriers occurs and the system simply wanders in its very high dimensional 
phase space through the few channels where no energy increase is required. We 
have already shown that at very low T, the system equilibration time, t(T), 
diverges exponentially. In that region, the typical observation time windows, 
tobs, are such that tobs/T <^ 1, and the system is in the early stage of its off 
equilibrium relaxation from its initial state. This is schematically the origin of 
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the flattening of Sa at very low T |T^. Notice that, in a system observed for an 
exponentially long time, i.e., for tobs/T ^ 1, the creep rate, Sa, would indeed go 
to zero. 

Interestingly, our model along with a saturation of the creep rate, S{T), also 
shows a saturation of the dissipation in the limit T — > 0. We show in Fig.^ the 
differential resistivity, p(T) = dV/dl, measured for the same value of the model 
parameters used in the calculation of the creep rate S{T) in Fig.|^ (the precise 
definition of V and / is postponed to the next section where we consider the I-V 
characteristics). The continuous curve superimposed to p(T) in Fig.^ corresponds 
to the linear fit p(T) = po + CpT. These results clearly show a saturation in p(T) 
at low T towards a finite value, in a way similar to the one recorded in S{T). 

The present scenario, where off equilibrium phenomena dominate the anoma- 



lous low T creep 
relaxation [I 



is supported by the experimental discovery of "aging" in the 



[57| , pT| , |5^ . In fact, strong discrepancies are found be- 

tween "quantum creep theory" predictions and the observed low T relaxation 
in many compounds |T^, 0, ^ . Interestingly, a unified picture begins to emerge 
of magnetic and transport properties. This will be discussed further in the next 
section. 
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Figure 9: The differential resistivity, p = dV/dl, in the ROM model is plotted as 
a function of the temperature, T, for N^xt = 10. The continuous superimposed 
curve is a linear fit. The saturation of p(T) for T — * well compares with the 
one of the creep rate, S{T). 
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4 The I-V characteristic 



Vortex flow in driven type II superconductors also shows strong memory and 
history dependent effects. Here, we outline the relations with magnetic proper- 
ties and propose a scenario for a broad set of these kind of phenomena ranging 
from "rejuvenation" and "stiffening" of the system response, to "memory" and 
"irreversibility" in I-V characteristics. In relation to recent experimental results 
|6l| , |6^ , we discuss in particular the nature of "memory" effects observed in the 



response of the system to an external drive, i.e., the I-V characteristic. Our 
model explains the peculiar form of such a "memory" of vortex flow at finite T 
and other "anomalous" properties such as the time dependence of critical cur- 
rents. The essential step is, again, to identify the relevant time scales in the 
dynamics. 

The system is zero field cooled and prepared by increasing Next at constant 
rate, 7, up to the working value (here. Next = 10). Then we monitor the system 
relaxation after applying a drive, I (due to an external current which induce a 
Lorentz force on vortices), in the ?/-direction. As in similar driven lattice gases 
p^ , the effect of the drive is simulated by introducing a bias in the Metropolis 
coupling of the system to the thermal bath: a particle can jump to a neighbouring 
site with a probability min{l, exp[— (A7^ — e/)/T]}. Here, ATi is the change in 
Ti. after the jump and e = ±1 for a particle trying to hop along or opposite to the 
direction of the drive and e = if orthogonal jumps occur. A drive / generates a 



voltage V p9|: 



Vit) = {vait)) (11) 

where Va{t) = v{t) is an average vortex "velocity" at time t ||15|. Here, v{t) = 
X J2iVi{t) is the instantaneous flow "velocity", Vi{t) = ±1 if the vortex i at time 
t moves along or opposite to the direction of the drive / and f j = otherwise. 

4.1 Memory effects in driven vortex flow 

We analyse a striking manifestation of "memory" observed in experiments where 
the drive is cyclically changed in the low T region A drive / is applied to 
the system and, after a time ti, abruptly changed to a new value Ji; flnally, after 
waiting a time t2, the previous / is restored and the system evolves for a further 



ts (see lower inset of Fig.|TO]). The measured V{t) is shown in the main panel of 
Fig.|l^ for T = 0.1. A flrst observation is that after the switch to Ji the system 
seems to abruptly reinitiate its relaxation approximately as if it has always been 
at Ji (see for example the dashed curve in Fig.|lO]), a phenomenon known as 
"rejuvenation" in thermal cycling of spin-glasses and other glassy systems 



The more surprising fact is, however, that for /i small enough (say /i -C /*, /* 
to be quantitatively deflned below) when the value I of the drive is restored the 
voltage relaxation seems to restart from where it was at ti, i.e., where it stopped 
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Figure 10: In the ROM model the voltage, V(t), is plotted as a function of time 
at T = 0.1 for a drive 7 = 1. As shown in the lower inset, after a time lag 
tl, the drive is abruptly changed to Ji for a time t2 and finally it is set back to 
its previous value. When / is switched to Ji the system seems to "rejuvenate": 
it suddenly restarts its relaxation along the path it would have had if / = Ji 
at all times (consider the continuous and dashed bold curves, corresponding to 
J = /i = 1 and / = Ii = 0.8, plotted for comparison). By restoring / after t2, 
the system shows a strong form of "memory" : if ^2 and Ji are small enough (see 
text) the relaxation of V{t) restarts where it was at ti. However, if t2 and Ii are 
too large, this is not the case, as shown in the upper inset. In this sense, the 
above is an "imperfect memory". 



before the switch to Ji (see Fig jTOl) . Actually, if one "cuts" the evolution during 
t2 and "glues" together those during ti and t^, an almost perfect matching is 
observed (see upper inset of Fig.|lO]). What is happening during t2 is that the 
system is trapped in some metastable states, but not completely frozen as shown 
by a small magnetic, as well as voltage, relaxation. These non trivial "memory" 
effects are experimentally found in vortex matter and glassy systems We 
call them a form of "imperfect memory", because they tend to disappear when 
the time spent at Ji becomes too long or, equivalently (as explained below) when, 
for a given t2, h becomes too high, as shown in the upper inset of Fig .110. 



4.2 History dependent I-V 

We now turn to the time dependent properties of the current-voltage character- 



istic. As in real experiments on vortex matter we let the system undergo a 
current step of hight Jq for a time to before starting to record the I-V by ramp- 
ing /, as sketched in the inset of Fig.|Tl|. FigJTl| shows (for T = 0.1) that the 
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Figure 11: The I-V obtained at T = 0.1 by ramping / after keeping the system 
in presence of a drive /q = 1 for a time to as shown in the inset. The response, 
V, is "aging" (i.e., depends on to) and, more specifically, stiffening: it is smaller 
the longer to. 



I-V depends on the waiting time to. The system response is "aging": the longer 
to the smaller the response, a phenomenon known as "stiffening" in glass form- 
ers ^ These effects are manifested in the violation of time translation 
invariance of two times correlation functions, already discussed. 

These simulations also reproduce the experimentally found time dependence 
of the critical current Usually, one defines an effective critical current, I^^-^ , 
as the point where V becomes larger than a given threshold (say Vthr = 10~^ in 
our case): one then finds that I^^^ is to and Iq dependent (like in experiments 
|6^ I^-^f is slowly increasing with to, see Fig.pl]). 

It is interesting to consider another current cycling experiment which outlines 



the concurrent presence of irreversibility and memory effects (see Fig.|T2|). The 
I-V is measured by ramping / up to some value Imax- Then I is ramped back 
to zero, but at a given value the system is let to evolve for a long time t^. 
Finally, I is ramped up again (see inset of Fig.|T2|). The resulting irreversible V{I) 
is shown in Fig.^. For I > 1^ the decreasing branch of the plot (empty circles) 
slightly deviates from the increasing one (filled circles), showing the appearance 
of irreversibility. This is even more apparent after t^: for I < the two paths 
are clearly different. Interestingly, upon increasing / again (filled triangles), V{I) 
doesn't match the first increasing branch, but the latest, the decreasing one: in 
this sense there is coexistence of memory and irreversibility. Also very interesting 
is that by repeating the cycle with a new 1^; (squares), the system approximately 
follows the same branches. This non-reversible behaviour is also found in other 



16 



glassy systems |Q. However, spin glasses, for instance, seem to show the presence 
of the so called chaos effects 0, The chaos effect is absent in our system as 
it is also in other ordinary glass formers 6i]. This kind of interplay between 
irreversibility and memory can be checked experimentally in superconductors and 
thereby assess the present scenario. 
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Figure 12: The I-V is measured at T = 0.1 during cycles of / (see also the inset): 
I is at first increased up to Imax (filled circles); along the descending branch of 
the cycle (empty circles), when I = 1^ (in the main panel the /^'s, for two cycles, 
are located by the arrows) the drive is kept fixed for a time = 10^ and then 
the cycle restarted; finally, / is ramped up again (filled triangles). For I > 1^, 
the first increasing ramp and the decreasing one (resp. filled and empty circles) 
do not completely match, showing irreversibility in the I-V. After waiting at 
Ijjj, a much larger separation is seen. However, by raising / again (filled triangles) 
a strong memory is observed: the system doesn't follow the first branch (filled 
circles), but the decreasing one (empty circles). Furthermore, in a cycle with a 
lower (squares), the same branches are found. 



4.3 Differential resistivity 

In Fig.|T3|, we plot the I-V recorded after ramping / at T = 1 (filled squares) and 
T = 0.1 (open circles). The low T I-V has the typical S shaped form experi- 
mentally found @, I, m, ig, irm, but, since we are above Tc, this is only an 
effect of short times of observation. The linear continuous functions in Fig.|13| 
are, in fact, the asymptotic I-V, i.e., those recorded after applying a drive / and 
measuring V in the long time regime (for t = 1.5 ■ 10^ in Fig.^). The same 
analysis applies to the differential resistivity, R = dV/dl, shown in the inset of 
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Fig.|T^. Also here one might think to see some characteristic regimes (defined, 



for instance, by the values J^, Ip of the inset of Fig.|T3D in the "short time" R{I)- 
They might be the off stationarity, finite temperature rests of crossovers between 
different plastic channels flow regimes typically found at T = 0, as discussed in 



1^, |T^, and references therein (see also |6^). Here, the linear behaviour 
of the asymptotic I-V indeed shows that the crossovers in the "short time" R{I) 
tend to slowly disappear with time, thus they cannot correspond to transitions 



among different driven stationary phases |T9|, [2^, ^ |61|, ^ . This conclu- 

sion holds despite the regular behaviour of Im and Ip with T also experimentally 
seen (for instance Ip seems to rapidly grow with T). An intrinsic structure in R 



can possibly be observed at sufficiently lower currents and temperatures ^8 



4.4 Voltage relaxation 

The natural step to understand the above observations is the identification of the 
characteristic time scales of the driven dynamics, which in the present model can 
be well accomplished. This we now discuss. Upon applying a small drive, J, the 
system response, V, relaxes following a pattern with two very different parts: at 
first a rapidly changing non-linear response is seen, later followed by a very slow 



decrease towards stationarity (see V{t) in Fig.|lJ for T = 1 and / G {1,2,3}). 
For instance, for / = 3 in a time interval At ~ 2 ■ 10^^, V leaps from about zero 
to AVi ~ 2 • 10~^, corresponding to a rate = |AV^/At| ~ 10~^. This is to be 
compared with the rate of the subsequent slow relaxation from, say, t = 2 ■ 10~^ 
to t = 10^, Tf ~ 10~^: Ti and Vf differ of 5 orders of magnitude. 

In agreement with experimental findings |^T], |66|, the slow relaxation of 



V{t) has a characteristic double step structure, which asymptotically can be well 
fitted by stretched exponentials |6^: V(t) oc exp(— t/ry)^. The above long time 



fit defines the characteristic asymptotic scale, ry, of relaxation. The exponent 
P and Ty are a function of J, T and N^xt (see inset Fig.|l^): in particular ry(/) 
decreases with I and seems to approach a finite plateau for I < I*, with /* ~ 0(1). 
In this sense, the presence of a drive / makes the approach to stationarity faster 
and has an effect similar to an increase in T. 

The outlined properties of Ty clearly explain the history dependent effects in 
the experiments previously considered. For instance, the "imperfect memory", 
discussed in Fig.|10|, is caused by the presence of a long, but finite, scale Ty in the 
problem: for a given Ii the system seems to be frozen whenever observed on times 
scales smaller than tv{Ii). Thus, if ^2 is short enough (t2 < ry(/i)) the system 
preserves a strong "memory" of its state at ti. The weakening of such a "memory" 



found for higher currents Ii in Fig. 10, is also a consequence of the strong decrease 



of Tv{I) with /. The phenomenon of "rejuvenation" (see Fig.|T^) is, in turn, a 
consequence of the presence of the extremely fast first part of relaxation found in 
V{t) upon applying a drive and of the above long term memory. The existence of 
the slow part in the V{t) relaxation also affects the "stiffening" of the response 
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Figure 13: The I-V is recorded by ramping / for the shown T. The continuous 
and dotted curves (resp. T = 1,0.1) are the asymptotic I-V, i.e., those where, 
for a given /, V is measured after waiting t = 1.5 ■ 10^ (see Fig.|T^). Inset The 
differential resistivity, R = dV/dl, for the same data of the main panel. The 
horizontal lines are from a linear fit to the asymptotic I-V. The characteristic 
values Im and Ip roughly locate crossover points in the "short time" R, which, 
however, disappear if t ^ oo. 



in the I-V of Fig.|lT], which is due to the non-stationarity of the vortex flow on 
scales smaller that ry. Actually, in Fig. 11, for a given I the value of V on the 
different curves corresponds to the system being probed at different stages of 
its non-stationary evolution. Finally, in brief, the fact that Ty(/) is smaller at 
high currents, /, and larger at small / (and T), is responsible for the surprisingly 
concomitant effects of irreversibility and memory of Fig.|12[ 

The origin of these time dependent properties of the driven flow, and in turn 
those of I-V's, traces back to the concurrent vortex creep and reorganisation 
of vortex domains. In fact, both with or without an external drive, the system 
evolves in presence of a Bean like profile (see inset of Fig.^) which in turn relaxes. 
An important discovery is that the characteristic times scales of voltage and 
magnetic relaxation are approximately proportional []T3|. This outlines that the 
non- stationary voltage relaxation is structurally related to the reorganisation of 
vortices during the creep (a fact confirmed by recent experiments 



5 Conclusions 

In conclusion, we showed that the replica mean field theory and Monte Carlo 
simulations of a schematic statistical mechanics lattice model |jl5[ for vortices in 
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type-II superconductors (a system of particles diffusing in a pinning landscape) 
offer a comprehensive framework of off equilibrium magnetic and transport prop- 
erties observed in vortex matter. Off equilibrium phenomena in many respect are 
known to show strong "universalities" f^, ^. In fact, here we considered either a 
mean field or a two dimensional version of the ROM model, which, interestingly, 
reproduces a very broad spectrum of experimental results. Molecular Dynamics 
simulations of more realistic systems, when existing, seem to confirm the present 
scenario |ff5|, 118], |19|, EOl, and, even if very demanding in the low T and high 



fields region, they can be an essential test for it. 
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Figure 14: The time evolution of the response function, V{t), for the shown values 
of the drive / (at T = 1 and Next = 10). In the asymptotic regime V{t) is well 
fitted with: V{t) oc exp[— (t/ry)^]. Inset The characteristic scale of relaxation, 
ry(J), as a function of /. For / — > 0, ry(/) seems to saturate to a finite value 
which implies Jc = 0. 

We have seen that the model shows a reentrant phase diagram in the field- 
temperature plane (B,T), analogous to what observed in vortex matter. More 
specifically, we discussed the off equilibrium, "aging", properties of magnetic 
creep. At low temperatures a crossover point is found, Tg{Nf.xt)-, where the sys- 
tem relaxation times become exponentially large. They seem to diverge a la 
Vogel-Tamman-Fulcher at a lower temperature, Tc{Nf,xt), where an "ideal" glassy 
transition point can be located. Magnetic creep changes its structure around 
Tg-. above Tg it shows power laws asymptotically followed by stretched expo- 
nential saturation; below Tg it is logarithmic. This corresponds to a change in 
microscopic vortex motion: from diffusive (above Tg) to strongly subdiffusive 
1l5[] . We showed that in the low temperature region the system is very far from 



equilibrium and its time correlation functions, no longer invariant under time 
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translations, have interesting dynamical scaling properties analogous to those of 
other "aging" systems. The above "off equilibrium" scenario also explains the 
surprising experimental discovery of a finite creep rate, Sa > 0, when T — > 
(previously interpreted in terms of "quantum tunnelling" of vortices [^) in our 
purely "classical" model. 

At not too high temperatures (but still well above Tg), magnetisation loops 
are typically found when M is plotted as a function of the applied field, including 
a definite "second peak" when the Ginzburg-Landau parameter is not too low. 
The "second peak" is associated with a new phase transition in the system. This 
can be difficult to see in experiments because samples can be significantly out 
of equilibrium, as shown by the dependences of the loops on the external field 
sweep rate. 

Vortex flow in driven type II superconductors also shows strong memory and 
history dependent effects. We have shown how creep and transport properties in 
driven media are related. We proposed a scenario for a broad set of these kind of 
phenomena ranging from "rejuvenation" and "stiffening" of the system response, 
"memory" and "irreversibility" in I-V characteristics, to history dependent criti- 
cal currents. 

The emerging unifying scenario of magnetic and transport properties in vortex 
physics has interesting relations with off equilibrium phenomena in other glass 
formers and complex fluids such as random magnets and supercooled liquids. 
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